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I.  Introduction:  Houlinear  Phanome am  In  ngctwgflgtte  theory 

Ths  subject  natter  of  this  survey  is  the  manifestation  of  nonlinear 
phenomena  in  a  wide  variety  of  problems  of  electromagnetic  theory.  The 
specific  problems  which  ve  intend  to  focus  on  include  shock  formation 
and  propagation  for  plane  electronagentlc  waves  propagating  into  an  isotropic 
nonlinear  dielectric  half -space,  the  possible  dissipative  effects  of 
anisotropy,  nonlinear  conduction,  dispersion,  and  relaxation  on  such 
wave  propagation,  and  the  development  and  propagation  of  singularities 
in  electromagnetic  wave  propagation  in  nonlinear  distributed  parameter 
transmission  lines. 

In  earlier  work  we  have  studied  the  problem  of  shock 
formation  in  non-dlssipative,  non-dispersive  nonlinear  dielectric  half¬ 
spaces  and  singularity  development  in  rigid,  infinite,  nonlinear  dielectric 
rods;  these  problems  have  been  treated  in  [1]  and  [2]  with  some  of  the 
results  described  in  SlI  below.  The  basic  achievement  in  [1]  was  to 
place  the  problem  of  shock  formation  within  the  context  of  the  general 
theory  of  homogeneous  quasillnear  systems  of  hyperbolic  conservation 
laws  and  thus  the  approach  differs  from  that  taken  in  earlier  work  by 
De  Martini,  Townes,  Gustafson,  and  Kelley  [3),  Broer  [4],  Jeffrey  [5], 
and  Korobeinikov  [6],  Kataev  [7]  and  Donato  and  Fusco  [8].  In  particular, 
setting  the  problem  within  the  context  of  the  theory  of  quasillnear  systems  . 
of  hyperbolic  conservation  laws  allows  for  a  relatively  simple  transition  to 
the  cases  where  the  media  is  either  anisotropic,  and  exhibits  a  smell  nonlinear 
conduction  current,  or  exhibits  dispersive  effects  or  relaxation  phenomena . 

It  is  shown  in  $11  that  in  each  of  these  cases  our  problem  may  again  be 
set  within  the  theory  of  hyperbolic  conservation  laws,  the  difference  with 
the  isotropic,  non-dispersive  case  being  that  inhomogeneous  terms  which 
are  possibly  dissipative  in  nature  make  their  appearance;  a  chief  goal  of 


the  work  described  In  Sllb  will  be  to  examine  conditions  under  which  shocks 
do  or  do  not  fora  and  in  this  regard  recent  work  of  Nishida  [9}  and  Slemrod 
[10]  on  dissipative  quasil inear  systems  and  of  Liu  [11],  Dafermos  [12], 
and  Dafermos  and  Hsiao  [13]  on  hyperbolic  conservation  laws  with  inhoaogenecity 
and  dissipation  nay  prove  useful.  Some  interesting  recent  numerical  studies 
by  Fisher  and  Bishel  [14]  at  Los  Alamos  on  the  Interplay  between  dispersion  and 
nonlinearity  for  Intense  plana  wave  laser  pulses  are  also  described  in  §XIb. 


In  Sill  we  consider  the  problem  of  shock  formation  and  propagation  in 

% 

distributed  parameter  nonlinear  transmission  lines  where  one  of  the  parameters 
C,  the  capacitance,  depends  on  the  voltage  drop  v  accross  the  capacitor; 
this  problem  has  been  considered  previously  by  Landauer  [15],  Jeffrey  [16], 

Riley  [42]  and  Kataev  [7].  However,  these  authors  considered  primarily  the 
formation  and  propagation  of  jump  discontinuities  in  the  first  derivatives  of 
current  and  voltage  in  the  transmission  line;  nore  were  these  problems  previously 
set  within  the  context  of  the  theory  of  quasil inear  systems  of  conservation 
laws  and  such  a  formulation  is  effected  in  Sill  where  it  is  conjectured  that 
for  sufficiently  small  Initial  gradients  of  current  and  voltage  shocks  will  not 
form  In  a  nonlinear  transmission  line  with  either  nonzero  linear  resistance  or 
leakage  conductance, whereas  such  shocks  will  form  if  these  Initial  gradients 
are  sufficiently  large.  Shock  formation  and  propagation  In  nonlinear  trans¬ 
mission  lines  is  of  some  considerable  Importance  for  as  Kataev  [7]  points  out 
experimental  tests  have  shown  that  such  electromagnetic  shock  waves  have  led  to 
technically  feasible  methods  of  producing  short  current  pulses  with  a  fast  rise 
time. 

Maxwell's  Equations,  in  the  absence  of  free  current,  free  charge,  and 
and  external  magnetization,  have  the  standard  form: 


-2- 


(1.1) 


«  -  curl  jg,  div  •  0 


■  curl  H,  dlv  D  ■  0 

V  #v. 


3 

in  sene  bounded,  open  domain  fl  «  R  ,  where  Is  the  magnetic  field,  JJ  the 

magnetic  intensity,  J |  the  electric  field,  and  Q  the  electric  Induction  field 

(or  electric  displacement).  The  latter  field  is  defined  by  $  =  +  £($)with  £(£) 

being  the  polarization  and  the  permittivity  of  free  space.  In  a  vacuum, 

£  ■  0,  H  where  Uq  Is  the  permeability  of  free  space  and  £q  >  0, 

_2 

Uq  >  0  satisfy  CqUq  ■  c  .  In  most  standard  texts  on  electromagnetic  theory 
only  linear  constitutive  relations  between  £  and  (and,  hence,  between 

and  £)  are  considered,  i.e.,  £(Jg)  •  (£-£q)£,  £  >  eq,  or  perhaps,  -  £  £ 

where  £  is  a  tensor  of  rank  four. 

In  [1912]  Voterra  propose  the  set  of  constitutive 


(1.2) 


£(x,t)  -  $  J£(x,t)  +  £(B(x,t)) 


$(x.t)  -  JiJ|(x,t)  +  jg(j$(x,t)) 

_oo 

which,  though  still  representing  an  a  priori  separation  of  electric  and  magnetic 
effects,  are  sufficiently  broad  enough  to  include  as  special  cases  the  linear, 
rigid  non-conductors  of  Maxwell  [1973],  the  Maxwell-Hopklnson  Dielectrics  [1897]: 

ft 

JJ(x,t)  -  £  J£(x,t)  +  J  $(t-T)J|(x,T)dT 


(1.3) 


-1. 


$(x,t)  -  p  $(x,t) 


£  >  0,  p  >  0,  $(t),  t  _>  0,  continuous  and  monotonlcally  decreasing,  and 

special  cases  of  the  holohedral  and  hemlhedral  Isotropic  dielectrics  Introduced 
by  Toupln  and  Rivlln  [21]  In  order  to  explain  the  phenomena  of  absorbtlon  and 
dispersion  of  electromagnetic  waves  in  dielectric  materials.  All  of  these 


Hr 
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special  cams  refer  to  linear  theories  and  ara  fully  discussed  in  the  recant 
monograph  by  this  author  [22]. 

As  is  clearly  indicated  in  Bloeabergen  [23].  nonlinear  propertlea  of  the 
conatltutive  relations 

(1.4)  ft  "  ft^ft*  ft  “  u($ft 

have  been  recognised  for  sone  tine.  i.e. ,  that  the  dielectric  constant  and 
magnetic  permeability  can  be  functions  of  the  field  strengths.  The  importance 
of  nonllnearly  in  problems  of  electromagnetic  theory  is  also  emphasised  in 
Kataev  [7]  who  points  out  that  "linearisation  of  the  parameters  of  a  media  is 
only  a  first  approximation.  For  example,  a  linear  relationship  between  the 
conduction  current  and  the  voltage  (Ohm's  law)  exists,  generally  speaking,  only 
if  a  number  of  conditions  are  observed:  constancy  of  temperature  and  composition 
of  the  medium,  and  absence  of  magnetic  fields". 

The  monograph  [23]  begins  with  the  constitutive  hypothesis  (1.4)  and  most 
recent  attempts  to  prove  focusing  of  polarised  electromagnetic  waves  propagating 
through  a  nonlinear  dielectric  media  have  been  based  on  constitutive  hypotheses 
of  this  form  with  y  -  const,  and  e(J£)  -  eQ  +  e2|  |j^|  |  ,  eQ  >  0,  e2  >  0. 

Similar  constitutive  relations  have  been  employed  In  the  work  on  shock  formation 
cited  In  [3]-[8],  e.g.,  in  [3]  the  nonlinearity  of  the  constitutive  relations 
is  reflected  In  the  consideration  of  a  material  with  an  intensity-dependent 
Index  of  refraction  and  it  is  shown  that  such  dependence  tends  to  distort  an 
optical  pulse  along  its  direction  of  propagation  thus  giving  rise  to  pulse 
self-steepening  (optical  shock  development) .  Other  Interesting  non-linear 
problems  for  dielectrics  which  have  been  treated  in  recent  years,  and  which 
are  based  on  constitutive  relations  of  the  form  (1.4),  include  the  work  of 
Kasakla  and  Venkataraman  [24]  who  use  the  method  of  characteristics  to  study 
the  asrly  phases  of  propagstlon  of  a  large  amplitude  electromagnetic  disturbance 
in  a  nonlinear  dielectric  slab  embedded  between  two  linear  dielectric  media. 


die  work  of  Venkataraman  and  Rivlin  [25]  who  present  a  method  for  calculating 
the  change  in  amplitudes  and  phases  of  harmonics  of  all  orders  for  a  plane 
electromagnetic  wave  propagating  in  a  nonlinear,  non-dissipatlve ,  isotropic 
dielectric,  and  the  work  of  Rogers.  Ceklrge.  and  Askar  [26]  who  apply  a 
modification  of  the  Bergman  integral  operator  method  to  a  hodograph  system 
describing  linearly  polarized  plane  electromagnetic  pulse  propagation  In  a  non¬ 
linear  dielectric,  thus  obtaining  closed  form  solutions  which  are  valid  up  to 
the  time  of  shock  formation. 

II.  Singularity  Development  In  Nonlinear  Materials 

a.  Shock  Formation  In  Non-dissipative,  Non-dlsperslve  Isotropic  Media 

We  assume  Maxwell's  equations  have  the  form  (1.1)  with  conduction  vector 
^  and  free  charge  density  p  •  0.  Our  constitutive  equations  have  the 
form  (1.4).  We  consider  a  linearly  polarized  plance  wave  of  the  form 
$  -  (0,  Ey(x,t),0),  Jg  -  (0,DY(x,t),0),  £  -  (0,  0,Bz(z,t))and  %  -  (O.O.H^x.t)) 

propagating  Into  the  region  x  >  0  occuppied  by  the  nonlinear  medium;  in  what 
follows  we  will  drop  the  subscripts  and  just  write  E  -  Ey,  D  *  Dy,  etc.  We 
also  set  e(0  -  e(0,C,0),  y(C)  -  li(0,0,c)»  VC  €  R1  and  assume  that 

(II. 1)  (Ce(C))'  >0,  (Cy(C))’  >  0*  at  least  VC  e  R1.  |c|  sufficiently 

small. 

(II. 1)  suffices  to  yield  local  hyperboliclty  of  the  resulting  quasllinear  system 
Substituting  the  assumed  fora  of  the  propagating  plane  wave  Into  (1. 1),  Maxwell's 
equations  reduct  to  the  quasllinear,  homogeneous,  hyperbolic  system 


with  a(E)  -  l/(e(E)E)',  b(H)  ■  1  /(y(H)H) ' .  In  all  of  the  previous  works 
which  dealt  with  this  problem,  i.e.,  [4]-[8],  the  system  of  evolution  equations 
has  been  expressed  in  terms  of  E  and  B;  the  system  (II. 2),  however  is  not 
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in  conservation  form.  To  rewrite  It  in  conaarvation  form  va  note  that 


(II.3)  D  -  fE  -J7FT*  B  -  fH  _* 

;0  Job«r 

so  that  ■  Et/a(E),  Bt  •  Ht/b(H).  Furtharmora,  In  vlav  of  (II. 1),  tha 
Integrals  in  (II. 3)  are  monotone  in  E  and  H,  respectively,  and  thus  in¬ 
vertible.  Therefore,  3E,  H  such  that  E  ■  E (D) ,  H  ■  H (B) .  The  system  (II. 2) 
then  assumes  the  form 


where  E(D)  -  X(D)D,  H(B)  -  Y(B)B  and  X(D)  -  l/fc(EO>)),  Yt»>  -  l#(H(B)). 
(II. 4)  has  the  form  of  a  hyperbolic  conservation  law  of  the  form 


For  nonmagnetic  materials  (v(J|I)  •  UQ) ,  H'( B)  »  1/uq  and  the  system  (II. 4) 

has  associated  real,  distinct  eigenvalues 


X.  -  -  *  -i;  /PW 

and  associated  characteristics  defined  by  the  equation  ^C***)* 

With  (II. 4)  we  associate  Initial  data  of  the  form  D(x,0)  ■  D  (x),  B(x,0)  - 

o 

Bq(x)  and  we  remark  that  a  priori  estimates  of  the  type  proven  by  Nlshida  [9] 
and  Slemrod  [10]  may  be  used  to  show  that  if  E'(C)  >  0,  VC  s.t.  |c|  is 
sufficiently  small  (a  consequence  of  the  definition  of  E  and  (II. 1))  and 
|su£  Dq(x)|,  Jsu£  Bq(x) |  are  sufficiently  small,  then  E’D(x,t))  >  0  for  as 

long  as  C1  solution  of  the  IV?  for  (II. 4)  exists.  Thus  the  characteristics. 


as  wall  as  tha  Riemann  Invariants  given  by 


(11.5) 


/E,(C)d8» 


4(D»B)  -  — ;  I  o  i^TcT  da. 


are  well-defined.  It  Is  a  simple  matter  to  show  that  ji ,  4  are  constant  along 
their  respective  characteristics  so  that  (11.4)  is  equivalent  to  the  diagonal 
system. 


(11.6) 


-  ^1(*»0|j  ■  *(x,(YtO  »t)  •  0 
4'(x, t)  -|£  +  Vx,t)lf  "  df  4(x2(5,t),t)  -0 

*1(Y,0)  -  Y.  *2^’°^  "  5*  e  R+ 


,  f  _ 

To  this  system  with  associated  Initial  data  A(x,0)  *  BQ(x)  +  — —  /£*(£)  ds 

etc.  we  may  apply  the  following  results  on  singularity  development 


(1)  [Lax,  [27].  If  D0(«)  Is  periodic  and  BQ  -  0  (where  we  extend 
from  R.+  to  R1  by  periodicity)  and  P'(0)  i4  0  (genuine  nonlinearity) 
then  finite- time  blow  up  must  occur  for 


*x(x.t) 


Bx(x,t)  + 


/E'  (D)  D 


-u0Dt(*.t)  +-r: 


x<X.t) 

Dx(x,t) 


i.e.  V^x  tjD  -  (Dx  Dfc)  must  blow-up  In  finite  time  and  a  shock  develops. 
Also,  if  the  Initial  values  differ  little  from  a  constant  4q,  4q,  then 

the  time  beyond  which  a  smooth  solution  cannot  be  continued  is  given  by 

(U.7)  *  ly'FW  /  (max  (x)|E"(0)|). 

R 

(11)  ([Klainerman  and  Majda  [29])  Suppose  that  Dq(*)>  Bq(*)  (and 

hence  4q(*)»  4q(*))  have  compact  support;  then  any 


solution  of 


K*  ■  4*“  •  0,  a.t.  *  C  oust  develop  singularities  In  finite  tine 

If  £'(£)  la  not  conatant  on  any  open  Interval. 

«e  now  describe  briefly  two  results  obtained  by  this  author  in  (2]  which  are 
based  on  the  above  theorems.  For  a  nonmagnetic  material  of  the  form 

D  -  <eQ  +  X0)B  +  X^2.  B  -  uqH 

we  have  E  -  XQD  +  X^2  +  0(D3)  3  E(D)  where  XQ,  Xj  are  related  to  the 

linear  optical  susceptibility  Xg  ^  the  £irst  nonlinear  optical 

2  3/2 

susceptibility  x1  by  Xo"Xl^£0+X0^  ’  *1  "  ”xl/  (eq+Xq)  *  As 

|E"(0)  |  -  2 1 Xj 1  +  0  we  may  take  BQ  -  0  and  EQ(0  periodic  and  apply  the 

results  of  Lax  [27].  In  the  MKS  system  eq,  xq»  Xt  a11  of  order  of 

1  -7 

magnitude  10  while  yQ  is  of  order  of  magnitude  10  .  Assuming  EQ 

9 

.  su^|e(x,0)  I  to  be  of  order  of  magnitude  10  volts/meter  (as  it  would 
be  in  an  intense  lightwave  or  laser  beam)  we  find  the  approximations 

t  «  VeQ +*& _ 1 

nax  X^/2  (ma£  Eq  )  (max  Eq  ) 

R  R 

vbeam  *  (y^i)  l^2(«ax|E0|)  1^2 

R 

where  s  is  the  maximum  distance  travelled  by  the  beam  into  the  media  until 
nax 

shock  formation  and  CQ  -  Uq^CEq+XqJ/x3  is  a  characteristic  constant  of  the 

particular  material.  The  above  expression  for  s _  compares  quite  favorably 

with  that  derived  in  [3]  especially  in  as  much  as  the  two  results  show  that 
8^^  la  inversely  proportional  to  the  largest  initial  gradient  of  the  beam 
(the  results  in  [3]  are  derived  in  an  entirely  different  manner).  One  of  the 
chief  benefits  of  writing  our  system  (II .4)  in  conservation  fora  is  that  we 
may  now  apply  standard  results  for  hyperbolic  conservation  laws  to  study  the 


formation  and  evolution  of  shocks.  The  Banklne-Sugoniomnt  conditions  Imply, 


in  particular,  the  possible  existence  of  two  shocks  with  speeds  s  , a.  given 
by 


(U.8)  vv  -  t  JL 

yo 

while  the  Lax  k-shock  conditions  [30}  (which  serve  as  an  admissibility 
criteria  to  single  out  a  unique  physically  realizable  weak  solution  of  IVP  for 
(11.4))  reduce  to  statements  that  for  s  ■  sr,  s  ■  either 


(II. 9) 


/ 


(e(E  )E  )’ 


(e(E+)E+) * 


or 


/  (e(E_)E_) 


>  >  - 


/  (e(E+)E+) 


where  E_,  E+  are  the  respective  values  of  E  just  in  back  of  and  in 
front  of  the  shock.  For  s  «  sr  it  is  shown  in  [1]  that  only  the  first 
of  the  Inequalities  is  applicable  while  for  s  •  only  the  second  Inequality 
is  applicable.  Furthermore,  it  is  shown  in  [1]  that  in  the  Important  special 
case 


e(E)  -  eQ  +  e2®2* 


£q  >0,  c2>  0, 


where  £"(©(5))  - 5 — ♦  E*  <I>(0 )  »  so  that  E"(0)  -  0  but 

(e0+3e2O 

E"(0(i;))  f  0  for  C  ^  0,  the  results  of  Klainerman  and  Majda  [29]  may 
be  applied  and  for 

s  -  s*  -  Uq1/2  (eQ+e2(E2  +  (E+  +  EJEJ)'1/2 
the  admissibility  criteria  imply  that 
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2  2  2  2 
2ET;  -  E  >  E,E  >  2E  -  E* 

2  2 

from  which  It  follows  that  E+  >  E_  for  the  shock  moving  to  the  right.  In 

an  analogous  fashion  we  show  in  [1]  that  the  admissibility  criteria  imply  that 
2  2 

E  >  E*  for  the  shock  moving  to  the  left.  As  the  local  electric  field  energy 
density  in  the  wave  given  by 

mao)  i  W  •  *  ■  K®2  + 


and  energy  must  be  dissipated  accross  the  shock , (Kataev  [7]), we  conclude  that 
the  shock  moving  to  the  right  with  speed  s  *  s^  is  not  physically  realizable. 

For  the  system  in  nonconservation  form,  i.e.*  (11.2),  it  can  be  shown 
(Kataev  [7],  Jeffrey  and  Korobeinikov  [6])  that  there  exist  simple  wave 
solutions  of  the  form 

(11. 11)  E  -  EQ(x  ±/a(D)b (H) t) ,  H  -  HQ(x  ±/a(E)b(H) t) . 

For  the  nonmagnetic  material  these  simple  waves  propagating  in  the  media 
assume  the  form 

(11.12)  E  -  EQ(x  t  p‘1/2^XE)t),  H  -  HQ(x  ±  y"1/2v£(E)  t) . 


It  was  essentially  deduced  by  Broer  [4]  by  a  direct  computation  based  on 
(II.  12j)  that  Ex(x,t)  »  as  t  -►  t*  -  2^  e*  if  3  e*  >  0  a.t. 


m 


VZ  e  ft1,  | ^  sufficiently  small;  this  condition  is  essentially 


equivalent  to  the  genuine  non-linearity  assumption.  Shock  development  in  the 
more  general  case(II.ll)  may  be  deduced  directly  from  the  work  of  Jeffrey 
and  Korobeinikov  [6],  i.e.,  it  is  possible  to  deduce  conditions  under  which 
E^,  Hx  blow  up  along  the  characteristics  defined  by  •  ±  /a(E)bdH) 

by  direct  differentiation  of  (II. 11)  and  use  of  the  fact  that 


(11.13) 


4(E,H) 

4(E,H) 
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and  Relaxation 


With  the  exception  of  some  rather  vague  remarks  In  Broer  [4],  concerning 
the  effects  of  dispersion  on  plane  wave  propagation  in  isotropic  nonlinear 
dielectrics,  the  numerical  studies  of  Fisher  and  Bishel  [14]  and  Shimizu  [38], 
on  the  interplay  between  nonlinearity  and  dispersion,  and  the  effort  made  in 
the  paper  of  DeMartini,  Townes,  Gustafson,  and  Kelley  [3]  to  build  an  exponential 
relaxation  process  into  the  light  pulse  equation,  most  treatments  of  wave 
propagation  in  nonlinear  media,  to  date,  have  ignored  the  effects  on  shock  forma¬ 
tion  and  propagation  of  anisotropy,  nonlinear  conduction,  dispersion,  and 
relaxation.  In  this  section  we  will  outline  an  exact  mathematical  formulation 
of  some  of  these  phenomena  emphasizing  chiefly  the  role  of  anisotropy  and 
nonlinear  conduction  and  indicate  how  recent  work  on  dissipative  and  inhomogeneous 
quasilinear  hyperbolic  systems  can  be  used  to  study  the  models  obtained. 
Electromagnetic  shock  waves  are  usually  quite  difficult  to  observe;  even  in 
the' case  of  ideal  isotropy,  with  no  conduction,  no  dispersion,  and  no  relaxation 
(instantaneous  response)  the  calculations  in  [1]  and  [3],  e.g.  the  expression 

for  s  in  §IIa,  indicate  that  whereas  shocks  always  form  the  physical 
max 

observation  of  these  shocks  depends  on  the  magnitude  of  the  largest  initial 
gradient  of  the  pulse.  Gradients  which  are  extremely  large,  but  still  within 
the  range  of  applicability  of  classical  physics,  are  needed  in  order  that 
shock  formation  occur  within  physically  observable  distances  of  several  meters 
(or  less).  In  [3]  the  authors  consider  an  initial  Gaussian  pulse  of  the 
form 

p(0,t)  -  pQ  exp(-4t2/t£) 

2  2 

where  p  ■  n^E  /8ir  is  the  approximate  energy  density  in  the  wave  and  nQ 


the  nonlinear  refractive  index;  t^  is  the  initial  width  of  the  Gaussian  in 
time.  The  expression  for  the  maximum  distance  travelled  until  shock  formation 
which  is  derived  in  [3]  is  of  the  form 


-1 

smax 


3v 

2Mt' 


min 


where  Vg  is  the  linear  velocity  of  propagation  of  the  wave.  Thus, 

varies  inversely  with  0^)^^  which  is  the  largest  negative  slope  of  the 

initial  pulse  in  time.  For  ^-switched  pulses  in  CS^  the  authors  [3]  Indicate 

that  t»  is  about  10  nsec  yielding  an  s  of  about  5  meters;  for  a  mode- 
4  max 

locked  laser,  however,  they  [3]  indicate  that  t^  can  be  less  than  10  ^  sec 
and  that  pulse  steepening  (shock  formation)  can  occur  over  propagation  path 
lengths  of  less  than  a  centimeter.  The  basic  conjecture  which  we  now  put 
forth  is  that  anisotropy,  nonlinear  conduction,  dispersion  and  relaxation  may 
act  as  dissipative  mechanisms  in  the  governing  evolution  equations.  That 
dispersion  and  relaxation  may  interfere  with  shock  formation  has  been  touched 
upon  in  [3],  [4],  and  [14]  via  appropriate  (numerical)  computations.  No 
study,  as  far  as  we  can  tell,  has  ever  touched  upon  the  influence  of  anisotropy 
or  nonlinear  conduction  on  shock  formation  in  nonlinear  electromagnetic 
materials,  even  though  no  material  is  perfectly  isotropic  and  certainly  no 
dielectric  is  a  perfect  nonconductor  (especially  when  interacting  with  an 

intense  light  wave  of  the  type  generated  by  a  laser) ;  we  therefore  begin  by 

\ 

formulating  a  mathematical  approach  to  the  influence  of  anisotropy  and 
conduction  on  plane  wave  propagation  in  nonlinear  dielectrics. 

We  again  consider  plane  waves  of  the  form  ^  -  (0,E(x,t) ,0) , 

H  •  (0,0,H(x,t))  and  assume  the  material  to  be  nonmagnetic  so  that  £  -pgg 
and  thus  £  «  (0,0,B, (x,t)) .  Our  assumption  of  anisotropy  takes  the  form 

£  and  are  not  parallel  in  the  media;  specifically 


l  “  t)),  i.e.. 


we  cake 


(11.15)  £(x,j|(x,t))  -  (600,  P(E(x,t)),0),  Se Cl(R+) 


with 


(11.16)  sup  1 6  (x)  |  <  6fl,  ,  sup.-r-<  <$.;  6Q  >  0,  6. 

xeR+  lx6R+  I 

Then  by  Che  definition  of  Che  eleccrlc  Indue Cion  field 


>  0 


(11.17)  £(x,c)  "  +  ^(x,J^(x.C)) 

-  (S(x),egE(x,t)  +  p(E(x,c)),0) 

-  (<5(x),  D(E(x,C)),0). 

Clearly  dlv  ^  ■  0  and  dlv  Q  -36/3x  so  ChaC  3<5/3x  represents  Che  free  charge 
density.  In  response  Co  Che  £  field,  and  by  virtue  of  the  existence  of  an 
effective  free  charge  density  36(x)/3x,  the  conduction  vector  £  can  be 
expected  to  be  nonzero  In  the  dielectric  with  £  "  cr(x,Jg)Jg  (nonlinear  Ohm's 
Law;  the  explicit  dependence  of  the  conductivity  a  on  position  is  consistent 
with  Inhomogeneity) .  In  view  of  our  constitutive  assumptions 

(11.18)  0,o(x,E(x, t))E(x,t)  ,0) 

where  a(*,0  3  a(» , (0,?,0)) ,  5  €  R1 


and  3/3t("|^x-)  +  V*>J  *0  so  that  the  equation  expressing  conservation  of 

charge  is  trivially  satisfied.  As  £  jf  J|j)  the  third  of  Maxwell's  equations 
must  be  modified  so  as  to  read 


i  +  3$/3t  “  7  *  J? 

while  the  fourth  equation  is,  of  course,  just  div  D  *  3<5/3x.  With  the 
assumption  of  plane  wave  propagation  these  equations  then  reduce  to  the  scalar 
quasllinear  system 


(11.19) 


S 


0(x,E)E 


.  as. 

3x 


as 

3t 
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Assuming  again  that  3D/3E  >0  so  that  3e  with  E(X,t)  -  £ (D(x,t)) 
nay  rewrite  (11.19)  in  conservation  font  as 


H  +  - -cf(x,E(D))E(l»  •  -  J(x.D) 

ft  +  I  ■  0 

where  £(x,D)  2  c(x,E(D))E(D) .  Clearly  (11.20)  Is  an  lnhonogeneous  conservation 
law  o£  the  form 


1  3B 


(11.20) 


3ui 

(11.21)  3t  +  3x  ■  gj/x.Jj) 


with 


D  (lM)\ 


The  system  (11.20)  may  be  compared  with  that  considered  by  Nlshida  [9],  and 
Slemrod  [10],  l.e.. 


(11.22) 


v  -  v  «*  0 
t  x 


vt-T(w)  ■  -dv 


(a  >  0) 


which  yields  the  damped  nonlinear  wave  equation 
*tt  *  ~t  ■  r(u)o- 

Initial  value  problems  for  (11.22)  have  been  shown  in  [9]  to  have  global 
C*  solutions  (wVv)  provided  the  Initial  gradients  wx(x,0),  vx(x,0) 
are  sufficiently  small  in  the  C*'  norm  while  shocks  have  been  shown  to 
develop.  In  spite  of  the  damping  factor  a,  if  these  gradients  are  sufficsntly 
large  polntwlse  [10].  Our  system  (11.20)  leads  to  the  wave  equation 

(11.23) 

As  £’(.,?)  -  E’(C)(a(’.E(c))-  E(C)$' (. ,E(;)) ,  and  E’(;)  >0,  y;  .  Rl 
(hyperbollclty) ,  we  con lecture  that  for  sufficiently  small  Initial  gradients 


Dtt  +  r(x,D)Dt  -  Eto)^. 
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& 


n»W*- 


D^(x,0) ,  E(x,0)  solutions  of  I VP  for  (11.20)  will  also  ba  globally  smooth 

if  £'(•,?)  >  0,  vC  €  I1<b>cf(«,y)  >  yo'  (*»Y)»  vY  e*l»  but  that  shocks  will 
davalop  if  thasa  initial  gradlants  ara,  pointwise,  suf f iclantly  large .  The 
analysis  in  [9],  [10]  is  basad  on  tha  equivalent  system  of  aquations  satisfied 
by  tha  Riemann  invariants  associated  with  (11.21).  For  our  system  (11.20)  tha 
Riemann  invariants  are  again  given  by  (11.5)  where  now,  for  the  sake  of 
simplicity,  we  normalize  and  take  ^  ■  1.  Then  (11.20)  is  easily  seen  to  be 
equivalent  to 

(11.24)  K'  -  -W  Z(x,D) ,  £(x,D) 


where  '  -  3/3t  +  X  3/3x,  '  -  3/3t  -  X3/3x,  X  -  WPW. 

As  -  |  /FTcT  dc,  a 6*  such  that 


D  -  S*(r-s)  (i.e.  j 


j  SFW  dC  -  B 
0 


B*(a)). 


Thus  (11.24)  assumes  the  form 


A'  ■  -p(x,*-4),  '  -  3/3t  +  X*(A-4)3/3x 

(11.25) 

4'  -  p(x,A-4),  '  -  3/3t  -  X*(A-4)3/3x 

where 

!p(x,A-4))  -  i /V  (0*(A-4)  •  L(x,6*(A-4)) 

Z(x,0*(/t-4))  -  a(x,E(0*(A-4)))E(3*(*-4)) 

X*  (A-4)  -  X(6*(A-4))  -  /E'(0*</t-4>>. 

Thus,  one  goal  of  further  work  in  this  area  should  be  to  investigate  the 
existence  or  nonexistence  of  global  smooth  solutions  for  I VP  for  the  system 
(11.25)  beginning  with  the  case  where  the  conductivity  is  homogeneous  so  the 
p  does  not  depend  explicitly  on  position.  In  this  regard  an  attempt  should 
be  made  to  generallza  and  extend  the  results  of  Nishida  [9]  and  Slemrod  CIO] 
for  systems  of  the  form 
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V  -  - 1  (A+d) .  a'-  - 

by  a»a«i<w<ng  the  behavior  of  A.,  4  along  their  respective  characteristics 
(in  our  problem  the  curves  defined  by  -  ±/E'  (6* (*-*))“  iX(8*(A-.A))  ■ 
X*(A-4).  Such  work  is  now  In  progress. 

He  note  here  that  we  could  equally  veil  work  with  E ,  H  Instead  of  D,  B. 

If  we  set  a(E)  -  3D/3E,  b(H)  -3B/3H  (we  do  not  now  begin  with  the  assumption 
that  the  media  is  nonmagnetic).  A(E)  -  l/a(E) ,  8(H)  ■  l/b(H),  then  (11.19) 


assumes  the  form 


(11.26) 


)+  /  o  A(E)\/e\  _  _  /  E(x,E)\ 

,t  (b(h)  *  o  /  y®  J  *x  \  0  / 


where  E(x,E)  -  -o(x,E)A(E)E 


The  characteristics  associated  with  (11.26)  are  given  by 


^  -  ±/A(E)8(H) 


and  the  left-eigenvectors  associated  with 


are  L  -  (1,±/A(E)/B(H)) 


Using  multiplication  by  the  componenets  of  the  first  left  eigenvector  to  combine 
the  equations  In  (11.26)  In  the  standard  fashion  we  find  that 


,  H'  _  _  (x.E)  _  _ 


vflC(E)  ^8(H) 


#(£)■ 


A(E)0(x,E)E  5  -A  (x  ,E) 


where  '  -3/3t  +  /A (E) 8(H)  3/3x  denotes  differentiation  along  the  characteristic 

4^.  -  /A(E)8(H).  Thus  if  we  define  the  Riemaniann  invariants 
dt 


!,H)  -  fE  v'StfST  dC  +  fH  « 

Jo  Jo 


'  (E.H)  -J 


W)  d; 


^b(t)  d; 


have 


K'  •  - A(x.E),  4"  -  -A(x,E) 
In  the  nonmagnetic  case 


(11.27) 


Jt 


and 


(11.28)  ft' •  -A(x,E)  ,  4K-  -A(x,E) 

with  '-3/3t  +  vXOO  3/3x,  '  -3/3t  -  3/3x. 

fE 

But,  by  (11.27),  ^  -  j  d?  so  3  8*  a.t.  E  ■  g*(-t+4). 

Setting  A(E)  •  /A(E) ,  we  then  have  as  a  consequence  of  (11.28) 


|*t  +  A*(*+4)*x  -  -f(x,*+4) 

4t  -  X*(4+4)Ax  -  -^(x.fcM) 

where 

!A*(A+4)  -  A(g*(A+4))  -  A(g*(*+4> 

¥(x,*+4)  -  +A(x,g*  (*♦*)) 

-  g* (A.+4)V'A(g*(/l+4))o(x , g* (-1+4) )  . 

Finally,  if  we  do  not  specialize  to  the  case  of  a  nonmagnetic  material  then 

¥  -  j‘  «•  ¥  ■  {]  « 

so  that  we  may  set 

E  -  g*(/t+A),  H  -  f*(4-4) 

A(E)  -  *£(17  -  /A  (g*  (*+*))  -  A*(/t+6) 

♦  (H)  -  -  4>*  (A.^6 ) 

and  in  place  of  (11.29)  we  obtain 


K  +  X*(*+4)$*(Jt-d)A  -  -?(x,*f4) 

(II. 30) 

■4t  ~  X*(*+aH*(4-A)  x  •  -4'(x.-t+A) 
with  ¥(x,*+*)  as  In  (11.29). 

We  could  also  approach  thasa  probisms  In  tarn  of  trying  to  apply  and/or  axtand 
results  on  tha  axlstanca  and  bahavlor  os  solutions  to  inhoaoganaous  hyparbollc 
conservation  laws  of  tha  fora  (11.21),  share  in  our  problam  ^(jj)  -  (g^)» 
when  u(H)  ■  const.  With  regard  to  such  problems  tha  existing  literature  is 
quite  small  but  soma  important  recant  work  (Liu  [11],  Liu  and  Li  [31],  Dafermos 
and  Hsiao  [13],  and  Dafermos  [12])  should  prove  to  be  of  Interest  for  our 
particular  problem.  Liu  [11],  in  fact,  constructs  weak  solution  of  (11.21) 
for  the  case  where  ^  -  ^(x,t)  is  an  n-vector,  £  is  a  smooth  n- vector¬ 
valued  function  of  jg  and  £  and  3,|/3^x  are  piecewise  continuous  n-vec tor- 
valued  functions  of  x  which  are  continuous  In  and  studies  their  asymptotic 

behavior  as  t  -*■  +».  It  is  assumed  that  /3jj  has  real  and  distinct 

eigenvalues  <  X2ty)  <  •••  <  Xq(jj)  for  each  jj.  The  analysis  in  [11] 

is  based  on  a  numerical  scheme  which  generalizes  the  Gliam  scheme  [32]  for 
hyperbolic  conservation  laws.  Among  the  results  proven  in  [11]  are  the 
following: 

(1)  when  the  X^ty(x))  i  -  l,2,...,n  are  nonzero  and  the  L1  norm  of 

g(x,jg(x>)  is  small  for  jg(x)  uniformly  close  to  the  initial  data  ^(x),  a 
global  solution  of  the  I7P  for  (11.21)  exists  and  tends  pointwise  to  a  steady 
state  solution  of  ■  £(x,;g). 

(li)  when  each  characteristic  field  is  either  genuinely  nonlinear  or 
linearly  degenerate  (Lax  [30])  the  solution  tends  uniformly  to  a  linear 
superposition  of  shock  waves,  refraction  waves,  travelling  waves  and  a  steady 
state  solution;  however,  these  waves  are  determined  by  the  values  of  the 
Initial  data  ^(x)  at  x  ■  f.  It  is  assumed,  in  the  general  case. 
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that  the  Initial  data  have  snail  total  variation. 

In  [31] ,  Liu  and  Li  study  I  VP  for  (11.21)  under  the  assumption  that  £(x.j£) 
has  coapact  support  in  x,  with  particular  emphasis  on  constructing  non- 
interacting  nave  patterns  for  conservation  laws  with  a  general  moving  source 
term  of  the  form 

(11.31)  +  -gg-  -  g(x-ct,ji). 

This  system  can  be  reduced  to  (II. 21)  under  the  change  of  variables 

"  c%  and  x  -*•  x  -  ct.  In  [13]  Dafermos  and  Hsiao  consider  I  VP 
for  systems  of  the  form  (11.21)  with  ^(x  j^(x,t))  -*•  ^(x,t,^(x,t))  and 
establish  the  existence  of  locally  defined  solutions  with  shock  waves. 

They  prove  the  existence  of  globally  defined  solutions  by  introducing  an 
appropriate  definition  of  dissipativeness  for  the  source  ^(x.t,^)  and  showing 
(via  a  combination  of  Gllmm's  scheme  for  conservation  laws  and  the  method  of 
fractional  steps)  that  its* effect  counterbalances  the  wave  amplification  due 
to  Inhomogeneity. 

Many  of  the  existing  results  for  systems  of  the  for  (11.21)  in  both  R*  and 

R11,  n  >  1,  with  an  emphasis  on  asymptotic  behavior,  are  sunarized  in 

Dafermos  [12]  while  problems  for  particular  forms  of  either  the  inhomogeneity 

j |  or  the  function  £(^)  have  been  treated  by  Liu  [33]  and  Ting  and  Wang  [34]. 

Concerning  the  possible  dissipative  effects  of  relaxation  on  the 

formation  and  evolution  of  shock  waves  in  nonlinear  electromagnetic  materials,  we 

should  be  Interested  in  replacing  constitutive  relations  of  the  form  (1.41) 

t 

by  assumptions  of  either  the  form  J£(x,t)  “  E(E(x,T))J£(x,t)  with 

-40 

t  ft  , 

(11.32)  E($(x,t))  -  eQ  +  I  e2(t-t)  |$r(x,T)  dx. 

and  |  -  (0,S(x,t),0),  H  -  (0,0,H(x,t)) ,  $  -  g 
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or 


£  -  (O.P(x.t).O)  with 


t  ft  2 

(11.33)  P(x,t)  -  P  (E(x,t))  -  *(t-T)E  (x,T)dT. 

J  ■  '00 

In  this  latter  case  ve  would  have  J^>  ■  (0,D(x,t),o)  end 

f  *  2 

(11.34)  D(x,t)  -  eQE(x.t)  +  j  \|Kt-T)E  (x.T)dr. 

It  then  follows  that  Maxwell's  equations  reduce  to  the  system  of  scalar 
quaslllnear  hyperbolic  Integrodlfferential  equations 


(11.35) 


|| +  ||  -  -l|/(0)E2(x.t)  -  J  i|/t(t-T)D2(x,T)dT 


3H  ,  3E 
3t  3x  -  0 


With  Riemann  Invariants 


X  -  E  +  — 1 —  H, 


and 


3/3t  + 


3/3x, 


4  -  E  -  — —  H, 

•V 


-  3/3t  - 


3/3x, 


we  obtain  the  following  system  of  nonlinear  integrodlfferential  equations  along 

\  dx  1 

X  the  characteristics  ■  ±  ■ —  . 


(11.36) 


X1  -  ~p|,(0)  Ot+4) 2  f  <|>t(t-T)  (*+4)2  (x,T)dT 

4'  -  -|^(0)  (*+4) 2  -J  j_  l|>t(t-  T)  (A+4)2  (x,T)dT. 


If  we  follow  the  constitutive  hypothesis  (11.32).  Instead  of  (11.33).  then 
we  obtain.  In  place  of  (11.38)  the  system 


-?n_ 


Ie0  *  |^2<t-T)E2(x*T)dTl  H  +  f  - 

(II. 37)  |  -e2(0)E3(x,t)  -  |jC  e2(t-T)E2(x,T)drj  E(x,t) 


2S  +  1I  -o 

3t  +  3x 


whose  characteristics  are  defined  by  the  equations 


dx 

dt 


+ 


i: 


e2(t-T)E2(x,T)dT]'1/2 


The  system  (11.36)  has  the  advantage  of  possessing  families  of  parallel  straight 
lines  as  characteristic  curves  In  the  x,t  plane  but  distortion  of  propagating 
waves  may  occur  due  to  the  nonlinear  source  terms  on  the  right  hand  sides  of 
these  equations;  an  analysis  paralleling  the  work  of  MacCamy  [35],  [36]  and 
Battori  [37]  on  materials  with  memory  may  be  possible.  The  form  of  the  con¬ 
stitutive  hypothesis  (11.32)  more  closely  follows  the  formulation  in  [3]  than 
does  (11.33)  and  a  numerical  treatment  in  [3]  indicates  that  a  rarefaction 
shock  may  still  build  up.  The  Importance  of  considering  relaxation  effects  in 
electromagnetic  media  la  emphasized  in  Kataev  [7]  who  indicates  that  equations 
(1.4)  "can  be  Invalid  for  material  media  (if  there  is)  a  lag  in  the  change  in 
induced  fields  which  occurs  in  response  to  a  rapid  change  in  the  field  intensities 
in  the  media.  This  lag  in  the  reactions  may  be  attributable  to  the  magnetic  or 
dielectric  viscosity  of  the  medium" 

Up  to  this  point  we  have  totally  Ignored  the  role  that  dispersion  might 
play  in  propagation  of  waves  in  nonlinear  electromagnetic  media.  In  a  linear 
dielectric  (£  -  £qX&)  where  the  dimensionless  coefficient  X  is  the 
polarisibility  of  the  medium  and  satisfies  the  wave  equation 


(11.38)  AE  -  e0P0(l+x)32E/3t2  -  0, 

so  that  the  wave  speed  is  given  by  c  ■  {£q1Jq(1+x)  l^2}  and  the  refractive 
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1/2 

Indue  n  ■  e/c  *  (1+X)  *  dispersion  and  relaxation  are  usually  introduced 

into  the  linear  theory  by  replacing  the  above  constitutive  relation  between 
£  and  by  the  Loren tz  equation 

<“•*»  7  r4  ♦  l  -  V*- 

For  a  wave  having  an  j|  field  with  monochromatic  frequency  u)  (11.41)  leads 

2  2  2  2 

to  the  well  known  result  that  n  -  1  +  X»q/  (u>Q  -  ut  ) .  One  way  of  introducing 
dispersion  into  the  nonlinear  optics  relations  (1.4)  for  an  isotropic  wedia 

(which  we  also  assume  to  be  nonmagnetic)  would  be  to  take 

« 

(11.40)  e  -  e(u>;E)  -  eQ(w)  +  e2(w)E2 

and  consider  the  propagation  of  plane  waves  in  which  E  has  the  form  (0,E,0) 
with, 

(11.41)  E(x,t)  ■  Re(Eg(x)e-iiWte*'CX) ,  k  the  wave  number. 

The  corresponding  reduction  of  Maxwell's  equations  is  a  straight-forward 

matter  and  will  not  be  pursued  here  (although  it  is  being  pursued  by  this  author 

within  the  scope  of  the  general  research  program  described  above). 

He  conclude  this  subsection  by  noting  that  both  the  numerical  studies  of 

Fisher  and  Bishel  [14]  as  well  as  the  earlier  numerical  studies  of  Shimizu 

[38]  have  shown  chat  shocks  can  form  on  the  leading  edge  of  pulses  in  a  Cs2 

laser  when  dispersion  and  nonlinearity  are  considered  simultaneously. 

Ill .  Electromagnetic  Shock  Wave  Formation  in  Nonlinear  Distributed  Parameter 
Transmission  Lines 

Shock  formation  and  propagation  in  nonlinear  transmission  lines  has  been 


studied  by  a  number  of  authors  ([7],  and  the  references  to  some  of  the  Russian 
literature  cited  therein,  [15],  [16],  [42],  and  [43]).  As  Jeffrey  [16]  points 
out,  the  method  of  analysis  that  is  used  for  linear  systems  does  not  easily 


generalize  to  the  study  of  nonlinear  distributed  parameter  transmission 
lines  since  the  customary  procedure  in  the  linear  theory  is  to  convert  the 
two  first  order  transmission  line  equations  into  a  single  second  order  equation, 
the  well-known  equation  of  telegraphy.  When  the  transmission  line  is  nonlinear, 
the  voltage  and  current  are  no  longer  solutions  to  the  linear  telegrapher's 
equation  and  superposition  of  solutions  no  longer  applies.  (Further  developments 
in  this  vein  require  that  the  nonlinearity  be  simplified  via  approximations  of 
various  kinds,  as  in  Ostrovskii  [43]). 

Landauer  [15]  and  Riley  [42]  both  exploited  the  formal  equivalence  between 
the  nonlinear  transmission  line  equations  and  the  one-dimensional  lsentropic  gas 
flow  equations  to  study  nonlinear  transmission  line  behavior.  In  fact,  Landauer' 
work  [15]  represents  a  direct  application  of  the  methods  of  gas  dynamics  [41] 
while  Riley's  [42]  contains  an  application  of  the  work  in  [15]  to  a  nonlinear 
transmission  line  Involving  a  voltage  dependence  capacitance,  but  Jeffrey  [16] 
notes  that  neither  of  the  author's  ([15],  [42])  "considered  the  problem  of  when 
a  continuous  wave  propagating  down  a  nonlinear  transmission  line  first  becomes 
discontinuous  and  forms  an  electromagnetic  shock  wave";  he  considers  in  [16] 
an  idealized  transmission  line  with  distributed  parameters  C(v)  the  capacitance 
R  the  resistance,  L  the  self -inductance,  and  G  the  leakage  conductance  per 
unit  length  of  the  line  (see  figure  below) 


Here  v(x,t)  is  the  voltage  at  a  point  Which  is  x  units  distant  from  an 
origin  taken  in  the  line;  we  also  denote  by  i(x,t)  the  current  at  (x,t) 
and  note  that  by  definition  C(v)  *  dQ/dv  where  Q(v)  is  the  voltage  dependent 
charge  per  unit  length  of  the  line.  Ohm's  and  Kirchoff's  laws  when  applied  to 
the  element  of  the  transmission  line  depicted  above  yield  the  system  of  equations 


(III.l) 


t  li  +  +  R1 

L  at  +  ax  +  Ri 


c(v)!?  +  i  +  Gv 


o 

o 


In  (III.l)  the  analysis  in  [16]  assumes  that  L,  R,  and  G  are  constant.  It 
is  clear  that  (III.l)  has  the  equivalent  form 


(III. 2) 


(VC  -('M 

\  v  / » t  \1/C(v)  0  /  \  v  /  ,x  I-  Gv  / 

\  C(v)  / 

of  an  Inhomogeneous  quasllinear  hyperbolic  system  if  C(v)  - 

the  associated  characteristics  are  clearly  defined  by  ^  1 

Ac(v) 


4^  >  o. 

dv 


independent  of  v,  C  * 


Cq  -  const.,  Jeffrey  concludes  [16]  that 


In  fact. 
For  C 


(III. 3a)  (i(~)1/2  +  v)  -  -(  —§3  i  +  ^  v) 

dt  co  u 


along  solutions  of 


dx 

dt 


(III. 3b)  <L  (i(-L_)  1/2  -  v)  -  -  (-5-  i  -  v) 

dt  co  D0 

along  solutions  of  4—  — 

dt  ^ 

and  shows  that  i,  v  both  propagate  along  the  linear  line  with  attenuation 
factor  -G/Cq  if  R/L  ■  G/Cq.  In  attempting,  however,  to  show  that  an 


-24- 


electromagnetic  shock  forms  in  the  line  with  C  ■  C(v)  >0,  K  ■  G  ■  0  he 


eqs.  (4.  la),  (4.1b)  claims  that  along  solutions  of 
respectively. 


That  A,  4  as  defined  by  (III. 4)  are  constant  along  the  respective  character¬ 
istics  defined  by  ^  ■  t  —  ^  —  follows  directly  from  (III. 2)  with  R  -  G  -  0; 

/LC(v) 

shock  formation  for  IVP  associated  with  (III. 2),  with  R  ■  G  *  0,  is  then  a 
consequence  of  the  results  of  Lax  [27]  for  the  genuinely  nonlinear  situation 
Q"(0)  0  or  Klainerman  and  Majda  [29]  for  the  case  of  compactly  supported 

initial  data  i(x,0),  v(x,0)  with  Q"(C)  JO  on  every  open  interval  of  Rl. 

Even  allowing  for  the  incorrect  form  for  the  Rlemann  invariants  in  [16],  shock 
formation  is  obtained  in  this  paper  only  for  waves  propagating  into  a  constant 
state  1q,  Vq  and  follows  from  the  fact  that  characteristics  of  one  family 
adjacent  to  the  constant  state  have  constant  slope  and  are  thus  straight  lines. 
Moreover,  the  shocks  obtained  are  weak  in  the  sense  that  the  fields  are  required 
to  be  everywhere  continuous  so  that"shock  formation"involves  the  development  of 
jump  discontinuities  in  the  first  derivatives  of  the  fields.  For  the  model 
(III. 2)  with  R^  +  0,  Jeffrey  shows  [16,  §5]  that  disconinuities  form 

if  either 

<«  Vf>0  <  0  ^ 
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1 


(ii) 


£_(-!r)  >  0  and 


3vy 


1  (-2£) 
x^  3v'o 


(»  +  -£) 

"  C0 


<  1 


where  the  zero  subscript  indicates  evaluation  at  the  constant  state  ahead  of 
the  wavefront  trace  and  the  superposed  ^  indicates  the  limiting  value  along  the 
wave  front  trace.  However,  as  Q  ■  Q(v)  and  C(v)  •  ^  >  0  the  expression 

fv 

Q(v)  •  C(C)d?  may  be  inverted  to  yield  v  -  t/(Q)  and  the  inhomogeneous 

v0 

system  (III. 2)  may  be  rewritten  in  conservation  form  as 


(  +  r  (,(,l),,  ‘  ~R1 

(III. 5)  < 

(  i? 

If  G  ■  0  (no  leakage  current  between  the  conductors)  then  by  direct  analogy 
with  the  damped  Quasilinear  system  (11.21)  [w  -*•  Q,  v  ■+■  i,  r  -*•  ^  v  and  a.  R] 

la 

it  should  follow  from  the  work  of  Nishida  [9]  and  Slemrod  [10]  that  C1 
solutions  of  I VP  for  (III. 5),  and  thus  also  for  (III.l)),  should  exist  globally 
if  ix(x,0)  and  (^(x.O)  are  sufficiently  small  in  thr  C*  norm. but  that 
shocks  should  develop  if  either  of  these  initial  gradients  is  sufficiently  large 
at  some  point  in  the  transmission  line;  this  contrasts  sharply  with  the  results 
of  Jeffrey  [16]  cited  above  where  discontinuities  are  proven  to  form  (albeit, 
weak  shocks)  if  ^xQ(y)g  <  0  and  the  wave  is  propagating  into  a  constant  state. 
Wbrk  is  now  in  progress  which  involves  analyzing  the  system  (III. 5),  with  G  f  0, 

by  re-writing  (III. 5)  in  terms  of  the  appropriate  Riemann  invariants  and  studying 
the  behavior  of  these  Riemann  invariants  along  their  respective  characteristics; 
in  this  sense  the  study  of  (III. 5)  is  expected  to  have  close  connections  with 


the  study  of  the  systems  (11.25),  (11.29),  and  (11.30). 

In  the  model  considered  in  [7],  Kataev  begins  with  the  equations 


(III. 6) 


iai  +  as.  o 

3t  3x 
3t  3x 


where  J  is  Che  leakage  current  and  4  the  magnetic  flux  per  unit  length 
between  the  conductors.  The  connection  with  (III.l)  is  that  4  -  4(1),  L(i)  ■ 
34/di  (the  differential  Inductance),  J  ■  J(l,v)  *»  G(i,v)«v  and,  as  previously 
assumed,  Q  -  Q(v);  in  (III.l),  where  R  i*  0,  it  has  been  assumed  that  L(i) 
and  the  conductance  G(v,i)  are  constant.  Kataev  [7]  then  specializes  to  the 
case  where  L(i)  -  Lq  »  const.,  C(v)  ■  Cq  -  const,  and  thus  obtains  from  (III. 6) 
the  system 


(III. 7) 


It  11  +  .  o 

L0  3t  +  3x  0 

co!ir  +  !£*  -c(liV),v 


dx 


whose  characteristics,  defined  by 


,  are  straight  lines  (we  assume 


that  LqC0  >0).  He  looks  for  solutions  of  (III. 7)  in  the  form  of  stationary  shocks, 
l.e.,  for  solutions  which  are  functions  of  w  ■  x  -  8t  where  6  is  the  velocity 
of  propagation;  in  this  case  1  ■  i(w),  v  ■  v(w)  and  (II. 7)  becomes 


(III. 8) 


v  -  vQ  -  8L0(i-i0) 


Iw  "  8C0  %  "  G(i(w),v(w))-v(w) 


after  integration  of  the  first  equation.  In  (III. 8),  1q,  Vq  refer  to  the 

initial  values  of  the  current  and  voltage.  Solutions  of  (III. 8)  are  then 
sought  for  which  i  ♦  1q  as  w  -+■  +°°  and  i  -*■  ig  +  ig  as  w  -*■  +*»,  where  ig 
is  the  amplitude  of  the  stationary  wave;  this  is  achieved  by  Integrating  (III. 8): 


_ dC  _ 

G[;,0Lo(;-io)  +  v0]*(L08(t-i)  +  vQ] 


and  then  imposing  the  boundary  conditions  so  as  to  conclude  that 
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(a)  for  iQ  <  i  <  iQ  +  1Q  and  v  -  0Lo(i-i.o)  +  vQ,  G(i,v)  >  0 

<b)  L  G(l!v)v  11  div*r**nt  for  v  "  ^0(i-i0)  +  vQ 
*0 

when  i  -  1q  and  i  ■  ig  +  £q. 

A  further  analysis  then  shows  that  (a),  (b)  require  that  as  we  approach  points 
on  the  profile  of  the  wave  G(ltv)  -*■  0  yielding  (in  Kataev's  notation)  at  each 
point  on  the  profile  a  relation  of  the  fora  v  a  y(i) .  Taking  this  last  relation 
into  account,  Kataev  then  rewrites  the  first  equation  in  the  system  (111.7)  in 
the  form 

(1IX.9)  |i  -  -<L0/^>  || 


for  which  (approximate ) simple  wave  solutions  of  the  form 


(III.  10)  i  *  I(x  -  *  t) 

L0  dl 

are  Immediate.  The  development  of  discontinuities  In  the  current  i  now  follows 
Immediately,  under  obvious  conditions  on  y(i),  via  a  direct  computation  of  the 
gradient  i^.  A  more  satisfactory  approach  to  the  study  for  (111.7) 
in  our  opinion,  consists  of  introducing  the  Riemann  invariants 


(III. 11)  A(i,v)  -  i  +  CQ  v,  4(i,v)  -  /L^i-/C^v 

which  satisfy  (along  their  respective  characteristics  defined  by  ^ 


the  equations 


(111.12)  < 

1‘ 


Gq(A+4,A-4)*  (A-4) 
^G0(A+6,A-4)»  (A-d) 


where,  Gn(A+A,A-4)  •  G(  —  (A.-K4) ,  —  (A-4 ) ) .  The  Inhomogeneous  terms  on  the 
°  2^0 

right  hand  side  of  (III. 12)  simplify  considerably  if  the  conductance  G(i,v) 
is  Independent  of  the  current  1. 
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